


X 



On the Erraticity in Random-Cascading a Model|^ 

ZHOU Yi-Fei LIU Qin TANG Ying CHENG ChunQ 
(Institute of Particle Physics, Huazhong Normal University Wuhan 430079) 



Abstract 

The erraticity in the random-cascading a model is revisited. It is found that in contrary 

^D ■ to the previous expectation, even in the pure single-a random-cascading model without putting 

CN ■ 

in any particle there exists erraticity behavior and the corresponding entropy indices do not 

vanish. This means that the dynamical fluctuations in a pure single-a model already fluctuate 

^S| , event-by-event. Models with multiple-a strengthen this fluctuation. Taking double-a model as 

example, the variation of the event-space fluctuation strength with the mixing ratio of the two 

> . 
I; — , a's is studied in some detail. The influence of particle number on the results when particles are 

cn 

CS| ' putted into the system is also investigated. 
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Since the observation of large local fluctuations in a high multiplicity nucleus-nucleus 
collision event induced by cosmic ray in lOSS^^^, similar phenomenon has also been ob- 
served in the accelerator experiments^^] . Since then, extensive attentions are paid to the 
non-linear phenomena in high energy collisions'^' . 

Due to the finiteness of particle number, the particle distribution fluctuates with re- 
spect to the probability distribution obtained from theoretical calculation. This is the so 
called statistical fluctuations. For a long time the popular idea was that the probability 
distribution obtained from theoretical calculation is smooth and the actual particle dis- 
tribution fluctuates statistically on top of the smooth probability distribution. However, 
is it really the case? Is it possible that there are non-linear dynamics which causes the 
dynamical probability distribution itself to be non-smooth? This question is still to be 
answered. Although a systematic quantum non-linear theory has not been constructed 
yet and a definite reply to this quesition is still impossible by now, non-linear phenomena 
have already been observed in high energy experiments'^ ^1 as stated above. 

In studying the non-linear dynamical fluctuations, besides the factorial moments aver- 
aged over event sample, the importance of the fluctuations of single-event moments in the 
event samples has also been noticed''^]. It is shown that these fluctuations are associated 
with the system's chaotic behavior. A parameter /i, called entropy index, is introduced to 
describe the chaotic behavior. An entropy index greater than zero is the signal of chaos. 
This method is called "erraticity analysis". 

The random-cascading a model is a simple model which can simulate the non-linear 
dynamical fluctuations. It has been used in Ref. [8-11] to study the dynamical fluctuations 
in event space, i.e. erraticity. Since the value of parameter a determines the strength of 
dynamical fluctuations, it was argued in Ref. [9] that there is no erraticity in the model with 
a single parameter a, and therefore a model with Gaussian-distributed a was proposed 
instead. 

In this letter the erraticity of a model is revisited. It is found that there is erraticity 
already in a pure single-a model. The corresponding entropy index is small but does 
not vanish. It increases with the increasing of a. The mixing of different a's makes the 
entropy index increases further. The dependence of entropy index on the mixing ratio 
is studied. It is shown through putting in particles that in all cases the entropy index 
calculated by particle number is bigger than that in pure a model and tends to the latter 
only when the particle number tends to infinity. This reflects the effect of statistical 
fluctuations. 



Let us first briefly recall the definition of erraticity and entropy indexl"^! . 

The single-event factorial moments F^^^ and probabilitic moments C^^^ are defined as 
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respectively, where a phase space region is divided into M equivalent sub-regions, pm and 
Um are correspondingly the probability and particle numbers in the mth sub-region. 

The fluctuation of moments F'^^ and C'^f' in the event space can be characterized by 
their event-space moments 

Ci? = kFf^'') I {Ft^Y . C^J = (C(^^n/(Cj^V- (3) 

If C^^'' or Cjj^g^ has power law behavior 

Cf,) oc M^'^"' (^) , ^(5) oc M^'"' (^) , (4) 



when M is big then we say that there is erraticity in the system. It can be characterized 
quantitatively by the entropy index 



Another method for calculating entropy index is to express Cp^q as 



(5) 



Cp,q = {<i>i^^n, (6) 

where 

$(e) = ^(e)/^^(e)^^ ^f^ = C(^V(Cf )• (7) 

Define 

S,= ($(^)ln<I.f)), (8) 

then in the scaling region, i.e. in the region where E^ depends on InM linearly, we can 
write 

In the following we will study the erraticity in the random-cascading a modelt^^'. In 
this model the M partition of the phase space region A is realized in u steps. At the first 



step, it is divided into two equal parts; at the second step, each part is further divided 
into two parts, . . . , and so on. The steps are repeated until M = AY/6y = 1^ . In this 
step-by-step partition how particles are distributed between the two parts of a given phase 
space cell is determined by the value of independent random variable iOyj^^ where jj,is the 
position of the sub-cell ((l<j;y<2'')), v is the number of steps. The value of the random 
variable lo^j^ is given by: 

1 + olt 1 — (xr _i 

^v,2'j~\ = , UJu,2j = , j = 1, •••, 2"- (10) 

where r is a random number distributed uniformly in the interval [-1,1], « is a positive 
number less than or equal to unity, which determines the value-region of the random 
variable u, and thus describes the strength of dynamical fluctuations in the model. After 
z/ steps, the probability in the mth sub-cell {m = 1, . . . , M) is 

Pm = ^lh ^2J2 ■ ■ ■ ^uj, ■ (11) 

Using these probabilities the probability moment C^^-* in each event is calculated according 
to (2), and the entropy index fig is obtained through Eq's.(2)-(9). 

Besides the single-a model described above, the double-a-mixing model will also be 
used later. It is constructed by letting riA/" events among the total A/" events have a = ai, 
and the remaining (1 — ri)A/' events have a = a2. 

In order to study the influence of limited particle number, we put a certain number 
A^ of particles into the M subcells according to the Bernauli distribution 

A^' 
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The particle numbers rim {f^ = 1, • • • , M) in every window are thus obtained and the 
event factorial moments F^'^^ can be calculated. The entropy index /i^ is then obtained 
through Eq's. (2)-(9). 

Eq.'s (5) and (9) are two equivalent methods for getting entropy index. We now discuss 
which one is more suitable for our purpose. 

What we are interested in is the case when the partition number M of phase space is 
very big. When the number of particles is small and that of sub-cells is big, it is probable 
that the particle number in every window is or 1. It can be seen from Eq.(l) that this 
will result in F^^^ = for all g > 1, and then lead to <I>^^'' and Ylq both vanish. Therefore, 
the second method fails in this case. So we choose the flrst method to get iiq. 



When using the first method, it can be seen from Eq.(5) that fiq is obtained through 
partial ly differentiating ■i/'^"^ with re spect to p near p = 1. Then how to choose the range 

of p value near p = 1 to fit the partial deriva- 
tive is also a question which we must answer. 
A too small range will lead to a wrong result 
due to the limitation of computer precision. To 
make the choise we show in Fig.l the ip-p plot. 
It can be seen from the figure that within the 
range 0.5 < p < 2 ip versus p is apprximately 
linear. Therefore, we take it as the right range 

for fitting the partial derivative dip /dp. 
Fig.l ip r^ p relation 
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In real calculation when a is very small (e.g. a=0.1) the value of single event prob- 
ability moment C'^^ is close to zero. This results in big errors in numerical calculation. 
In order to solve this problem, we multiply the value of single event probability moments 
by a large number, e.g. 1000. This extra factor will appear in both the denominator and 
the numerator in the formula for calculating Cpg, Eq. (3), and will not effect the final 
results. 

The calculating results of single a model show that there exists erraticity. As example, 
in Fig. 2 is shown the variation of Cp^q with M for q;=0.3. It shows a typical characteristic 
of erraticity. The entropy indices /i2 in this model for various values of a are plotted at 
the right side of Fig.3. It can be seen that /i2 > 0, indicating the existence of erraticity. 
The value of entropy index /i2 increases with the increasing of a. 





Fig. 2 The variation of Cp^q ~ M 
for a = 0.3 



Fig.3 Entropy index for single 
a model 



On the left of Fig. 3 is shown the variation of entropy index with the number of 
particles when particles are putted into the model. It can be seen that when the number 
of particles is small the corresponding entropy index is much bigger than that of the pure 
a model without particle. The entropy index decreases with the increasing of the number 
of particles, tends to the value of the pure a model when the number of particles tends 
to infinity. 

The reason why erraticity exists also in the single-a model is that, in the expression 
Eq. (10) for the elementary probability in the model, there is a random number r, causing 
the window probability to fluctuate inspite of the fact that the values of a are the same 
in different events. 

It is clear that if the a's are different in different events, the fluctuation will be strength- 
ened. We take the double-a model as example to study this situation. 

The two a's are taken to be a\ = 0.2, 0:2 = 0.8. There are in total A/" events anioung 
which riJ\f events have ai = 0.2 and the remaining (1 — ri)A/' events have ai = 0.8. 

In Fig. 4 are shown the logio/U versus the number of particles n for ri = 0, 0.1, 0.3, 0.5, 0.7, 
0.9, 0.95, 0.98, 1, repectively. It can be seen that the entropy index fi firstly increases with 
the increasing of ri starting from ri = 0, but when ri is big (near to unity) it turns to 

decreasing. The reason is that ri = and ri = 
1 correspond to single-a model with a = 0.2 
and 0.8 respectively. The entropy indices for 
these two cases are in general smaller than that 
of the cases with the two a values mixed. The 
only exception is when only a very little frac- 
tion ( < 0.1) of big a (0.8) is mixed in. In this 
case, the double-a-mixing model will have a 

n smaller entropy index than the single-a model 

Fig. 4 Entropy index versus . , _ n 8 

multiplicity for different 

double-a model 

The variation of entropy index in double-a-mixing model with the mixing ratio ri is 

plotted in Fig. 5 for ai, 0:2 = 0.1, 0.9; 0.2, 0.8; 0.3, 0.7; 0.4, 0.6, respectively. The increasing 

of fi with Ti when the latter is not very big and the decreasing of fi when ri near to unity 

can be seen in all the cases. 

In order to show the fine structure near ri = 1 the curves for ri,r2 = 0.1,0.9 (case-A) 
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and 0.4,0.6 (case-B) are magnified and shown in the small figure inside Fig.5. The crossing 



(A) 



„(B) 



„(A) 



„(B) 




of these two curves is understandable. This is because al < a\ while 0:2 > « 

Therefore, when ri is small (1-ri is big) /i*-^-* > 
/i'-^-', but when ri is big (1-ri is small) /i*-^^ > 

In summary, it is found in this letter that 
even in the pure single-a random-cascading model 
without putting in any particle there exists er- 
raticity behavior and the corresponding entropy 
indices do not vanish. This observation is in 
contrary to the previous expectation. It means 
that the dynamical fiuctuations in a pure single- 
a model already fiuctuate event-by-event. 
This fiuctuation is strengthened when different values of a are mixed in the model. 
Taking double-a-mixing model as example, the variation of the entropy index with the 
mixing ratio of the two a's is studied. It is found that the entropy index fi firstly increases 
with the increasing of ri and then decreases. The maximum is located near ri = 1 when 
ai < a2- 

The infiuence of particle number on the results when particles are putted into the 
system is also investigated. In all cases the entropy index calculated by particle number 
is bigger than that in pure a model and tends to the latter only when the particle number 
tends to infinity. 



Fig.5 Entropy index in double-a 
model as function of mixing ratio 
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